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A combined experimental and detailed numerical srudv was conducted on the interaction between 
chemically inert solid particles and strained, atmospheric methane/air and propane/air laminar dames, 
both premixed and non-premixed. Experimentally, the opposed jet configuration was used with the addition 
of a particle seeder capable oi operating in conditions ot varying sravitv. The particie seeding system was 
calibrated under both normal and micro gravity, and a noticeable grautation.il effect was observed. Fiarne 
extinction experiments were conducted at normal gravity bv seeding inert particles at various number 
densities and sizes into the reacting gas phase. Experimental data were taken tor 20 and 37 u nickel alloy 
and 25 and 60 u aluminum o>ade particles. The experiments were simulated bv solving along the stagnation 
streamline the conservation equations of mass, momentum, energy, and species conservation for both 
phases, with detailed descriptions of chemical kinetics, molecular transport, and thermal radiation. The 
experimental data were compared with numerical simulations, and insight was provided into the elfects 
on extinction of the fuel tvpe. equivalence ratio, flame configuration, strain rate, particle type, particle size, 
particle mass delivery rate, the orientation of particle injection with respect to the flame, and gravity. It 
w'as found that for the same injected solid mass, larger particles can result in more effective flame cooling 
compared to smaller particles, despite the fact that equivalent masses of the larger particles have smaller 
total surface area to volume ratio. This counterintuitive finding resulted from the fact that the heat exchange 
between the two phases is controlled bv the synergistic effect of the total contact area and the temperature 
difference between the two phases. Results also demonstrate that meaningful scaling of interactions be- 
tween the two phases mav not be possible due to the complexity of the couplings between the dynamic 
and thermal parameters of the problem. 


Introduction 

The physical mechanisms of two-phase reacting 
flows are poorly understood compared to those ot 
pure gaseous flows. While substantial amount of 
work has been devoted to sprays, less attention has 
been given to the details of reacting dusty flows. 
Even chemically inert particles can alter the flam- 
mability and extinction limits of combustible 
mixtures through strong dynamic and thermal cou- 
plings. The presence of Stokes drag, phoretic, and 
gravitational forces on the particles affect the dy- 
namic interaction between the phases. Also, as the 
solid phase possesses larger thermal inertia relative 
to the gas. substantial temperature differences be- 
tween the two phases can develop. 

The stagnation flow configuration has been used 
in the past for the study of spray and particle flows, 
as it offers advantages both in experimentation and 
modeling [1-5]. Numerically, the need for a hybrid 
Eulenan-Lagrangian approach has been identified 
bv Continillo and Sirignano [1] and has allowed for 
the prediction of the phenomenon of droplet flow 


reversal [2.3], which has been observed experimen- 
tally [3]. Gomez and Rosner [4] have conducted a 
detailed studv on the particle response in opposed- 
jet configurations to determine thermophoretic dit- 
fusivities. Sung et al. [5] have numerically studied 
dustv flows in opposed-jet configurations, with em- 
phasis on the dynamic coupling. An extensive list ot 
pertinent studies can be founcl in Ref. [6]. In pre- 
vious studies on dusty flows [4,5], the thermal cou- 
pling between the two phases as well as the effects 
of particle mass delivery rate and gravity were not 
considered. In Ref. [6] such effects were partially 
addressed for premixed flames through detailed nu- 
merical simulations of the two-phase stagnation flow 
configuration. 

The objective of the present work was to conduct 
a combined experimental and detailed numerical 
simulation in order to assess the role of inert parti- 
cles on flame extinction. The results were used to 
derive physical insight into the mechanisms that con- 
trol the dynamic and thermal interactions between 
the two phases. Gravity was included as a parameter, 
as its effects can be quite complex. 
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Fic. L : i) Genera] experimental configuration, (b) Schematic of the particle seeder. 


Experimental Approach 

The experimental configuration, shown in Fig. la, 
includes the use of two counterfiowing jets exiting 
from two opposing 22 mm diameter burners sepa- 
rated by a 14 mm distance. Fig lb depicts a sche- 
matic of the particle seeder, which utilizes a piston 
feeder as proposed in Ref. (7). Attached beneath the 
bottom burner is a piston that feeds the particles into 
the Bow at a constant rate. Chemically inert alumi- 
num oxide and nickel alloy particles were used 
tthese dames are too cool for nickel ignition). Par- 
ticle sizes were determined to be within 2-5 u of the 
nominal size for 90% of the particles. The particle 
mass delivery is determined by both the piston speed 
and the flow rate. The gas flow enters the top of the 
piston shaft through sixteen 1 mm diameter holes 
equally spaced around the shaft, which locally in- 
crease the gas velocity and improve the entrainment 
of particles into the flow. This design allows for seed- 
ing under both normal (1 g) and micro gravity (ag). 
However, it should be obvious that the particle 
pickup is strongly affected bv gravitational forces. 
The particle seeder was calibrated bv seeding the 
particles into the air flow lor a specified time and 
their subsequent collection in ultralight bags at the 
burner exit. The delivered mass of the particles was 
determined by measuring the mass of the bag before 
and after the experiment bv using a high-precision 
scale. Calibration in pg was performed on board 
NASA’s KC-135 plane. 


Pre mixed and non-premixed flame extinction ex- 
periments were conducted in 1 g by vaning the par- 
ticle type, size, and mass delivers' rate as well as the 
gas phase equivalence ratio, fuel tvpe, flame config- 
uration, and strain rate. In ail these expen ments, the 
particles were seeded from the lower burner. 

Two types of premixed flame configurations were 
studied. The first included the use of the symmetric 
twin-flame that results bv impinging on each other 
two fuel/air jets of identical composition. For this 
configuration, the particles preferentially cool the 
first flame that they encounter and mav or may not 
have a chance to direcdv cool the second flame, de- 
pending both on the ability of the particles to pen- 
etrate the gas-phase stagnation plane iGSP) and on 
the particle’s thermal state as they reach the second 
flame. The second configuration used is a single 
flame that results by impinging a fuelyair jet on an 
opposing air jet. Given tnat the particles are seeded 
from the lower burner, the extinction of single flames 
was studied by stabilizing them on both sides of the 
GSP. If the single flame is stabilized below the GSP, 
then it is directly cooled by the particles. If the single 
flame is stabilized above the GSP. only large particles 
would possess the inertia to cross the "stagnation 
plane and directly affect the flame. 

In the non-premixed flame experiments, a mixture 
of fuel and nitrogen gas was supplied from the bot- 
tom burner, while the air jet was supplied from the 
upper burner. 

For all cases, the flames would be first established 
at conditions close to the extinction state. Then, the 
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FlC. 2. Particle seeder calibration cun-es lor 25.0 u ulu- 
ininum oxide particles in normal anti micro gravity. 

local strain rate, based on the maximum ve loci t\' gra- 
dient in the hydrodynamic zone, was measured by 
using laser Doppler ve loci met rv «. LDVV Subse- 
quently, the piston seeder was turned on. ieedingthe 
particles at a constant rate. The fuel How rate was 
then decreased very slowly, until the flames were 
extinguished due to the svnergistic etfect ol the 
strain rate and the heat loss to the particles. In order 
to isolate the effect of particle size on flame cooling, 
the piston speeds were set to deliver equal amounts 
of mass of different size particles The accuracy of 
the data is limited by the accuracy of calibration of 
the sonic nozzles, which is 0.5%. and the accuracy 
of LDV, which is 1 cm/s. The accuracy ot calibrating 
the particle mass flow rate was found to be less than 
10 %. 


Numerical Approach 

A set of equations along the svstem centerline has 
been developed for both phases [6]. A quasi one- 
dimensional set of equations was derived for the gas 
phase along the stagnation streamline, similar to 
those in Ref [S], The particle equations were for- 
mulated for each particle independently, valid for 
small particle number densities. 

For inert particles, the gas phase continuity and 
species equations are identical to those of Ref [8], 
The gas-phase momentum equation, however, was 
modified bv adding a term representing the force 
exerted bv the particles on the gas. The gas-phase 
energy equation was also modified by including 
terms describing the conductive/convective/ radia- 
tive heat exchange between the two phases. The ax- 
ial particle momentum equation includes the Stokes 
drag, thermophore tic, and gravitational forces: the 
svstem configuration is assumed to be vertical so that 
gravity acts in the axial direction. The particle energy 


equation includes the contributions o i conductive/ 
convective heat exchange between the two phases 
and the radiative energy exchange between the par- 
ticles and the surrounding gas. Finally, a conserva- 
tion equation is used to determine the particle num- 
ber density. 

The solution was obtained bv simultaneously in- 
tegrating the entire system of equations for both 
phases. For the integration of the gas-phase equa- 
tions, a Eulerian frame of reference was used, while 
the particle equations were integrated in a Lagran- 
gian frame of reference in order to properlv describe 
particle reversal (1,6). Detailed kinetics (9) were 
used, and the code was integrated with the CHEM- 
KIN ( 10 1 and TRANSPORT [11] subroutine pack- 
ages. 

Special care had to be taken in regions of particle 
reversal, as particle number density is singular (al- 
though integrable) at the points of reversal (6). In 
the vicinity of particle stagnation plane sPSPV the 
particle velocities are small so that at equal time 
steps, the Lagrangian spatial resolution is substan- 
tially greater than the Eulerian one. Thus, conver- 
gence mav not be possible if singular v alues ot par- 
ticle number densities are mapped from the 
Lagrangian grid onto a more coarsely resolved Eu- 
lerian grid, causing discontinuities in p rope rty values 
and thus in the phase interaction terms of the gas- 
phase momentum and energy equations. It was 
found that this problem can be resolved if the Eu- 
lerian grid is finely resolved around PSPs. 


Results and Discussion 

Experimental Results 

Flame extinction experiments were conducted on 
the extinction of methane and propane flames tor 
the following four cases: 

1. PremixecL symmetric across the CSP twin- flames 

2. PremixedL single-flame stabilized below the CSP 

3. PremixecL single- flame stabilized above the GSP 
-4. Non-premixed flames stabilized below the GSP 

In order to determine the particle mass delivery, 
the seeder had to be calibrated for all particle types 
and sizes. An example of the calibrations is shown in 
Fig. 2 for 25 p aluminur(0>xide particles. It can be 
seen that the particle mass flow rate increases line- 
arly with both gas flow rate and with the piston 
speed. Comparisons between the calibration curves 
at 1 g and ug reveal that for a given piston speed and 
gas flowrate, the particle mass flowrate at pg is con- 
sistently higher as no gravitational forces oppose the 
particle motion. For the same reasons, calibration 
below 150 cm 3 /s flow rate was not possible in I g, 
although easih’ achieved in pg. 
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(b) Equivalence Ratio at Extinction (d) 

^FlC. 3. Experimental data on twin premised flame extinction. Injection velocities at nozzle exits: fa) 46 4 cm/s :bi 
00.4 cm/s. (c) and (d) 33.5 cm/s. 


Flame extinction experiments performed in 1 ^ 
can be seen in Fig 3, which depicts comparisons of 
the cooling effects of 20 and 37 p nickel alloy and 
25 and 60 aluminum oxide particles on lean pre- 
mixed propane/air and methane/air twin flames. All 
flames shown in each frame experience the same 
strain rate. Note that in all cases, the data for the 
larger particles fall below those for the smaller par- 
ticles, indicating that for the same particle mass de- 
livery rate, larger particles extinguish stronger 
flames. This result is counterintuitive, as for the 
same particle mass the smaller particles expose more 
surface area to the gas phase, which should allow 
more efficient heat transfer. 

The results of Fig. 3 also reveal that for the same 
conditions, the extinction equivalence ratios are 
higher for propane flames, as their Le > 1 so that 
their extinction is more susceptible to strain rate and 
heat losses. It was also found that for the same con- 
ditions, the nickel alloy particles extinguish stronger 
flames compared to aluminum oxide. This may be 
explained by a combination of two effects. First of 


all, nickel alloy particles have about twice the specific 
heat of aluminum oxide. Second, the nickel allov 
particles are nearly spherical, while the aluminum 
oxiderSre quite angular As a result, the aluminum 
oxide particles have a larger surface- area- to- volume 
ratio, much like smaller spherical particles. 

Figure 4a depicts a comparison of flame extinction 
results obtained for the case of single methane/air 
premtxed flames stabilized below the CSP and twin 
methane/ air premixed flames. Comparisons of single 
methane/ air flames stabilized below and above CSP 
are shown in Fig. 4b. Recall that particles are in- 
jected from the bottom burner in all cases. Notice- 
ably, stronger flames are extinguished in the bottom 
flame configuration compared to the twin flame case 
with almost the same strain rate and particle mass 
delivery rate. This is physically reasonable, as twin 
flames are adiabatic, while single flames suffer from 
heat losses to air stream. The results of Fig. 4b in- 
dicate that the particles are able to extinguish 
stronger flames in the bottom flame configuration 
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Fie. 4. Experimental data on extinction of (a) twin and single flames stabilized above CSP and ib) single flames 
stabilized below and above CSP Injection velocities at nozzle exits 53.4 cm/s. 


compared to the top flame configuration. In the bot- 
tom (single) flame configuration, particles reach the 
flame region with temperatures nearly equal to the 
ambient temperature at which thev were injected. 
However, in the top (single) flame configuration, 
particles cross the stagnation plane and travel 
through the region of hot combustion products, the 
downstream of the top flame. As they move against 
an opposing flow, their velocities axe also reduced 
and their residence times within the hot product in- 
creases. Thus, their temperatures are raised close to 
the flame temperature dv the time they reach the 
reaction zone of the top flame. As a result, the par- 
ticles cannot effectively cool the top flame, and only 
relatively weak flames can be extinguished. 

Figure 5 depicts results of non-premixed flame ex- 
tinction. Fuel/nitrogen mixtures were supplied from 
the bottom burner and air was supplied from the top 
burner. For these reactant concentrations, the 
flames were stabilized below the GSP. that is. on the 
particle injection side. Again, stronger flames are ex- 
tinguished by the larger particles in all cases. 

Numerical Results 

Numerical simulations were conducted by inject- 
ing 25 and 60 /i aluminum oxide particles into twin 


flame methane/air flames. It was observed that, as 
in the experiments, 60 a particles cool the flames 
more effectively than the 25 u particles. 

In order to understand this behavior, the detailed 
flame structure was analyzed. Figure 6a depicts the 
spatial variation of the gas and particle phase tem- 
peratures obtained for both cases; in this and all fol- 
lowing figures, the particles are injected from the let t 
boundary, which corresponds to the lower burner. 

Note that there are two mechanisms that affect 
the heat transfer between the two phases. The first 
is the total particle surface area available for heat 
transfer, ana the second is the temperature differ- 
ence between the two phases. For the same mass of 
particles delivered, small particles possess larger to- 
tal surface area. Thus, from this respect, the total 
area mechanism favors the small particles. Howev er, 
because of their thermal inertia, large particles de- 
velop larger temperature differences with the gas 
phase. Thus, the temperature difference mechanism 
favors the large particles. It is the competition be- 
tween these two factors that determines which size 
cools more effectively. 

Figure 6b depicts the spatial variation of the heat 
transfer between the two phases for the two cases 
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Frc. 5. Experimental data on non-premixed flame extinction, with i 


i injection velocities at nozzle exits 33.5 cm/s. 


being considered. It can be seen that 25 ft particles 
remove more heat from the flames in the preheat 
region, and their temperature closelv follows the gas 
phase (Fig. 6a). However, the 60 u particles stav rela- 
tively cool until they reach deeper into the flame, 
resulting in increased heat loss within the reaction 
zone. Note that in the reaction zone of the bottom 
flame, the temperature difference between the two 
phases for the 60 ft case is of the order of the ratio 
of the total surface areas of the two sizes. Thus, the 
temperature difference factor dominates the area 
factor within the reaction zone, demonstrating why 
large particles can extinguish stronger flames. 

Figure 7 depicts additional numerical results 
which reveal that the physics can be even more com- 
plex Here, the variation of the maximum flame tem- 
perature with the strain rate is shown for methane/ 
air premixed flames at an equivalence ratio of 0.643 
seeded by 20 and 37 u nickel alloy particles. In both 
cases, the delivered particle mass density is 0.065 
mg/cm . Note that the large particles yield larger 
reductions in the gas-phase temperature up to a cer- 
tain strain rate (here about 140 s’ 1 ). Bevond this 
strain rate value, smaller particles exhibit a better 


cooling efficiency. (It should be noted that ali of the 
experiments presented above were below this criti- 
cal strain rate.) This is easily understood as. at high 
strain rates, the large inertia of the large particles 
rapidly transports them through the flame, leasing 
little time to remove much heat from the flames. On 
the other hand, small particles have less inertia as 
they approach the stagnation plane and reside longer 
within and thus can effectively cool the flame 
A quite complicated effect of both gravity and the 
particle mass delivery rate was foundrFig. 8 depicts 
the spatial s-ariation of the gas and particles velocities 
for a methane/air flame with equivalence ratio of 
0.584 seeded by 60 aluminum oxide particles, at 0 
g and 1 g, and for two different mass delivery rates, 
that is, 0.5S and 34.44 mg/ s. In the 0 g case, particles 
penetrate deeper into the upper jet since gravity 
does not reduce their velocities as thev approach the 
GSP. Fig. Sa and b depict the cases with 0.58 mg/s 
mass delivery rate. In both cases, particles undergo 
reversals and finally reach equilibrium at the staV 
nation plane for the 0 g case (Fig. 8a) and at a lower 
elevation, where Stokes drag balances the grarita- 
honal force, for the 1 g case (Fig. 8b). Fig. Sc and d 
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Fie. 6. Numerically determined la) spatial variation of the gas and particle phase temperatures and <b) interphusial 
heat transfer, for twin, premixed methane/air flames with an equivalence ratio of 0.6. seeded with 25 and 60 u aluminum 
oxide particles. 



FlC. 7. Variation of numerically determined maximum 
gas-phase temperature with strain rate for twin, premixed 
methane/air flames with an equivalence ratio of 0.645. 
seeded with 20 and 37 u nickel-allov particles. Particle 
mass density (n p *>n ? . where m p is a single particle massi 
delivered at the exit is 0.065 mg/cirr*. 


depict both results for mass delivery rate of 34.44 
mg/s. In the 0 2 case, however, panicles behave 
much as in the lower mass delivery rate case. Notice, 
however, the significant mass deliverv rate etfect be- 
tween the two 1 g cases, as shown in Fig. Sb and d. 
In the first case tFig. 3b), the mass delivery rate is 
so small that has minimal effect on the particle tra- 
jectories. But in Fig. Sd. the mass deliverv' rate is 
large enough to remove significant amounts of en- 
ergy from the lower flame. Thus, the thermal expan- 
sion and the resulting gas velocities are so small that 
Stokes drag cannot overcome gravity, and the parti- 
cles fall back into the bottom nozzle. 

These results demonstrate that the effects of par- 
ticle injection are very complex, so complex as to 
suggest that any simple scaling of the problem muv 
be all but impossible. For example, one might be 
tempted to suggest a scaling centering on a Stokes 
number based on the injection velocity iperhaps 
with the addition of a Froude number to capture the 
gravity effects). However, such a scaling misses a 
mvriad of effects due to the strong coupling between 
the particle velocities and the thermal effects. For 
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example, as it is based on the injection velocities, 
such a scaling does not account tor the presence or 
absence of particle reversals that are strongly af- 
fected bv the gas accelerations in the flame expan- 
sion zones. This can be clearlv seen bv comparing 
Fig. Sb and d, which show that the smaller gas ac- 
celeration near a weaker flame can cause particles to 
fall back into the burner l Fig. Sd) rather than un- 
dergo a series of reversals (Fig. Sb). (Also, for the 
same reasons, the gravity effect on the location of 
the reversal would not be captured bv a Froude 
number) To describe these eflects, one would also 
have to include thermal parameters, reflecting the 
strong effect of cooling on the Arrhenius behavior 
of the chemical kinetics i which cannot be scaled due 
to its nonlinear dependence on the temperature). 
One would be tempted to scale the particle mass 
delivery rate {apparent in Fig. S) with the particle 
volume, but such a scaling would not reflect the par- 
ticle size effects apparent in Figs. 3-7, as it is pos- 
sible to have the same solid-phase concentration for 


a variety of sizes. For that matter, there does not 
appear to be a convenient length scale with which 
the diameter could be scaled that would properlv 
reflect the effects of diameter on flame cooling. Con- 
sequently, no dimensional scaling of the problem 
makes itself apparent, and if possible would include 
a very large number of parameters that it could be 
scarcely said to simplify the problem. 

Concluding Remarks 

The effect of inert particle clouds on the extinction 
of strained premixed and non-prem ixed flames was 
studied experimentally and numerically in the coun- 
terflow configuration. The effects of fuel ripe, equiv- 
alence ratio, flame configuration, strain rate, particle 
ripe, size, delivery' rate, injection orientation, and 
gravity were independently considered. 

The experimental results revealed the surprising 
observation that large particles can cause more ef- 
fective flame cooling compared to smaller particles. 
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The numerical simulations reproduced this behavior 
and showed that while similar masses of small par- 
ticles possess more total surface area, the large par- 
ticles can establish large temperature differences 
with the gas phase within the reaction zone and. 
thus, more effectively cool the Hame. At high strain 
rates, however, the opposite trend was observed as 
large particles are quicklv transported through the 
Hame. The effect ot gravity was found to be quite 
involved, coupling the strain rate, the Hame tem- 
perature. and the particle mass delivery rate. 

Finally, the complex coupling between the various 
parameters is so complex that meaningful scaling 
capturing all pertinent phvsics is not possible. Thus, 
caution is recommended in using traditional scaling 
arguments to describe phenomena related to react- 
ing dusty Hows. 
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